Temperature dependent conductances of deformable molecular devices 
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Transport through a molecular device coupled to a vibrational mode is studied. By mapping it 
to the Yu- Anderson model in the large contact broadening limit, the zero bias electric and heat 
conductances are evaluated non-perturbatively. These exhibit a step from their T — value to 
half of its value as T increases due to the opening of the inelastic scattering channel. The spectral 
function exhibits the Franck-Condon suppressed multiphonon steps. The Wiedemann-Franz law is 
satisfied at low and high temperatures, but is violated in between. Relations to experiments are 
discussed. 
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Motivation. Molecular electronic devices are based 
on electron transport through individual molecules, and 
have been proposed for many years as possible candi- 
dates for future nanoelectronic circuit elements [l], Q. 
They triggered intense research due to the fundamen- 
tal challenge of quantum transport in nanoscale sys- 
tems and by the possibility of tailoring molecular scale 
structures [3, \M- O nc °f their most appealing features 
is the importance of local Coulomb interaction between 
quasiparticles and scattering off local vibrational modes, 
which are more pronounced due to the less effective 
screening mechanisms at small length scales. 

A wide range of characterization methods have been 
developed for collecting information about the micro- 
scopic properties of the molecular junctions and probing 
its excited states, including conductance histogram tech- 
niques B| ; Coulomb blockade @, Q , conductance fluctu- 
ation Is I, shot noise Q or superconducting subgap struc- 
ture [lOj measurements, operating mainly at low tem- 
peratures. An especially important information about 
molecular junctions is the study of vibrational modes 
with point-contact or inelastic electron tunneling spec- 
troscopy [E [lH, where the low-temperature zero-bias 
conductance of the device is perturbed by the excitation 
of molecular vibrations with an appropriate bias voltage. 
Indeed, experiments on single molecular devices includ- 
ing fullerenes^, 0], carbon nanotubes[I3|. benzene [l3j 
and simple molecules (H 2 , CO, H 2 0) connecting metal- 
lic electrodes [E EH, have revealed the influence of vibra- 
tional degrees of freedom supplemented by theoretical 
simulations of the vibrational modes 0, EH, EE EE 17 . 



IE HUES 

Whereas the above low-temperature techniques pro- 
vide a large amount of information about the ground 
state properties of the device and small perturbations 
due to the excitations of the inelastic degrees of freedom, 
any room temperature application requires the knowl- 
edge of the temperature sensitivity of the device and 
the impact of the strongly enhanced inelastic scatter- 



ing processes should be considered Therefore, in this 
work, we focus on the temperature dependent trans- 
port through molecular devices by studying the electrical 
and thermal conductance of a model system with non- 
perturbative analysis, valid for arbitrary electron- vibron 
couplings. 

The model and its mapping. To model a molecular 
transport junction, we consider the molecular Hamilto- 
nian as 
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H mo i = g v Qd + d+ — 
Zm 



-Q 2 



(1) 



describing a single electron level coupled to an Einstein 
phonon. The level is tunnel coupled to leads as 

H a = J2 <k)c+c k + J2 Vk (c+d + d+c k ) + e d+d, (2) 



where Vk is the hybridization parameter. The contact 
broadening is given by A = irpo(V k 2 ) fs- Let us con- 
sider symmetric contacts with non-interacting quasipar- 
ticles, in which case it suffices to consider a single effective 
contact [2 1| with energy dispersion e(fc) = vpk and den- 
sity of states (DOS) po = 1/2ttvf- First, we diagonalizc 
Hq, and express the molecular Hamiltonian in terms of 



its eigenfunctions. This is achieved by introducing 22 1 



and 
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where L is the size of the sample. The explicit form of 
r]k,k' is irrelevant for our discussion [^j. From now on, we 
restrict our attention to the cases when the broadening 
by the contacts is such 1J| that the DOS of the contact 
and the device are slowly varying on the scale of the 
phonon energy ujq and temperature T . This is realized in 
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the case A » (T,ljq) and A > |eo|, when the junction is 
represented by a resonant level model, thus Vk\f~L w T> = 
y / 2v F A/(e 2 . + A 2 ) is fc independent, which are the basic 
conditions for the mapping to hold. These are often sat- 
isfied under realistic conditions[l4ll5l|. Other cases were 
studied in Ref.dll Hit [Til . In terms of the transformation 
in Eq. (|3|), the total Hamiltonian Hq + H mo i is rewritten 
as 



which holds true even in the presence of vibrational 
modes, i.e. the zero temperature zero frequency den- 
sity of states remains unchanged. This occurs because 
at T = 0, the incoming electrons experience a frozen 
Fermi see and no bosons in the oscillator, thus no phase 
space for scattering. Identical results are obtained for the 
Kondo model. This leads to the conductance at T = 
as 



H = ^e(fc)a+a fc +.gQ*+(0)*(0)- 
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(5) 



where g = v 2 g v and ^(x = 0) = J2 k ajJ_\/~L. Eq. 
([5l) is known as the Yu- Anderson model[23[, describ- 
ing conduction electrons interacting with a local bosonic 
mode. The model is solved by bosonizing the fermionic 
field [H 0] as = exp[iv^$(x)]/v / 2lra, and the 

resulting effective model of one dimensional coupled har- 
monic oscillators [261 reads as 
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— Q 2 , 



(6) 

where g is the phase shift caused by g, $>(x) stems from 
the bosonic representation of the fermion field. The vi- 
brational mode softens, and the damped frequencies are 
given by 24 1 



G(T = 0) 



h A 2 



(11) 



which ranges from perfect transmission (G(0) = e 2 //i) 
for | eo | <C A to a decent suppression of the conductance 
to ~ 0.7 - 0.8e 2 /h for |e | < A. Eq. (JTTJ is expected to 
hold for low transparency junctions, beyond the validity 
of our mapping as well. The heat conductance satisfies 



T^O T 



L G{T = 0), 



(12) 



where Lq = (7rfcs/e) 2 /3 is the Lorentz number, thus 
the Wiedemann-Franz law is fulfilled. At high tempera- 
tures (T >Tnax(\ii}p±\,\<Jp+(jJp—\/T), but still obeying to 
T« A, the T-matrix for the 4" field reaches its univer- 
sal value, T = l/i2irpo. This halves the corresponding 
Green's function at x = in this high T region as 



G\j,(uj) — Gvj, + G%TG% 



iirpo 



(13) 



lo p ± = — iT ± 



; 2 (i-r/r 2 )-r 2 , 



(7) 



where T 2 = ttoj 2 /4W <C ojq W, W is the bandwidth of 
the conduction electrons, and T = ■n(gp$) 2 /2m for small 
g, and the model becomes unstable for T > T-2- The ex- 
plicit dependence of T on g v cannot be determined by the 
bosonization approach[24j. The T < Ti w T 2 (l — T^/lOq) 
region corresponds to undcrdamped phonons, while the 
overdamped response shows up at Ti < T < T 2 with two 
distinct dampings. 

Conductances at low and high temperatures. The elec- 
tric (G) and heat conductance (k) through the molecular 
transport junction in the wideband limit are given by [21 



G 



h I duj 



Tto 2 



ImG d (w), 



(8) 



where / is the Fermi function, and Gd{u) is Green's func- 
tion of the electron on the molecule. From Eq. ([3]), 



G d (w) = 9 2 G*(lu) 



(9) 



at x = 0. This basic relation allows us to determine 
the properties of localized electron from the Yu- Anderson 
model. At T = 0, one can derive a Fermi liq uid relation 
for the Green's function of the $ field [Til. 27 1 at x = as 



G$(w = 0) = —iirpo, 



(10) 



with G^ = —iirpo, which determines the conductance at 
high temperatures as 



G(T » w ) 



G(T = 0) 



(14) 



The phonon state is populated by many bosons at high 
temperatures, and every incoming electron scatters off 
them inelastically with increasing probability even in the 
weak coupling limit. The inelastic scattering cross sec- 
tion reaches its maximal value [27| at this temperature 
range for arbitrary electron-vibration coupling. The elec- 
trons dephase completely, and can scatter forward and 
backward with equal probability (=1/2). Consequently, 
the T-matrix takes the above universal value, and the 
conductance halves. Similar conclusions about halving 
the spectral weight were reached for scattering on a dy- 
namical boundary condition [2jj|. 

In the above high T limit, the heat conductance also 
satisfies k(T > w )/T = f^G{T > u ), and the 
Wiedemann-Franz law holds. Therefore, both the elec- 
tric and heat conductances are expected to exhibit a step 
from their T = value to half of its value with increas- 
ing temperatures, and the Wiedemann-Franz law holds 
in the two limits. In between the two limits, G and K fol- 
lows a different temperature dependence: the heat con- 
ductance drops faster with the temperature, and their 
ratio stays below the universal Lorentz number, hence 
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FIG. 1: (Color online) The electric and heat conductance 
(blue solid and red dashed line) through a molecular junction 
and their ratio (black dashed-dotted line) is shown as a func- 
tion of the temperature for V — 0.8r2, W = IOluq uj±/uq — 
±0.44 - i0.06. The high (T >Recj p± ) and low (T <Rew p± ) 
temperature parts are universal. The inset shows the evo- 
lution of the spectral function. p\z(uj) = —lmG\s/(uj)/n for 
T/lo = 0, 0.2, 0.4, 0.6 and 1.2 from top to bottom. 



the Wiedcmann-Franz law is violated, as is shown in 
Fig. [TJ The minimum of their ratio occurs roughly at 
3T =Recjp±, which can be used as a rule of thumb to es- 
timate the value of the renormalized phonon frequency. 
We mention, that by further increasing the temperature 
to leave the T <C A regime, Eq. JT]) practically decou- 
ples from the conduction electrons. In this regime, the 
conductances exhibit a second step to zero. 

Green's function from bosonization. To study the 
crossover between the high and zero temperature limits, 
we use the exact result for the Yu-Andcrson model ob- 
tained from bosonization 24| 23] valid for arbitrary tem- 
peratures, whose derivation is sketched below. From Eq. 
©, the electron Green's function on the molecule is re- 
lated to that of the ^ field. The real time dependence 
of the latter at the impurity site is obtained in the Yu- 
Andcrson model as 



G*(f) 



; Q(t) 



£ lim <{*( 7 M),* + (7'y,0)}) 



4 L — ' x,y— >0+ 
7,7'=± 



8tt 



Y, (exp[C 7 , r (t)]+exp[C y , 7 (-t)]), (15) 

7,7'=± 



where C 7 . 7 '(i) = lim Xty _ >0 + 47r($(7x, t)4>( 7 '?/, 0) - 
[<P{jx,t) 2 + $(7'?/,0) 2 ]/2) - ln(a), and a ~ l/W is the 
short distance cutoff in the bosonized theory, 7 and 7' 
denotes the sign of the spatial coordinates x and y. The 
expectation value C 7i7 ' (t) at bosonic Matsubara frequen- 
cies can be evaluated from the bosonized Hamiltonian, 



Eq. © as 

C 7i7 ' (w m ) 



4|w. 



1 l r (sgn(o; m ) + 7 )(sgn(cj ro ) - 7') 

2 (\u} rn \+lLU. 



ID. 



P-) ' 

(16) 



where u) m = 2rmrT. The first term is responsible for the 
1/t decay of the local fcrmionic propagator, while the 
second one stems from the interaction of electrons with 
the oscillator. For 7 = 7', this correction term vanishes. 
C 7i7 '(i) is evaluated by analytically continuing the Mat- 
subara frequencies to real ones in Eq. (|16|). which defines 
the spectral intensity, and then following Ref. 0, we 
arrive to the desired correlator. By defining the integral, 



A(t) 



dhJ- 



-i cxp(— iu>t) 



4T 



1 — exp(— uj/T) (u> — LUp+)(uj — Up-) ' 



(17) 



which follows the derivation of the position autocorrela- 
tor of a harmonic oscillator coupled to a heat bath[2g], 
and can be evaluated in closed form using the hyperge- 
ometric functions, the local retarded Green's function at 
finite temperatures follows as 

Gv(cj,T) = -i*^(l + I ph (w,T) + exp(A(0))), (18) 



where 



dt Te lf"!h mexp(A(t)), (19) 
smn(7TJ t ) 



whose T = limit was analyzed in Ref. 1271 giving 



ImG* = —iirpQ at T = u> = 0. In the T 3> |w p+ c<jp_|/r 
limit, I p h vanishes, and the imaginary part of the prop- 
agator is ImG^ = — i7rpo/2, which is the desired result. 
Plugging Eq. Q18|) to j9]), the conductances are evaluated 
from Eq. © , which are shown in Fig. [1] for the full T 
range, which agree nicely with the analysis of limiting 
cases. The inset shows the local spectral function, the 
Franck-Condon steps[l8|,[2(| arising from multivibron ex- 
citations are smooth due to the significant phonon damp- 
ing. Notice its explicit temperature dependence through 
the phonon occupation number. The power of the exact 
solution manifests itself in comparison with lowest order 
perturbation theory (LOPT) [13, , which neglects mul- 
tiphonon contributions (becoming dominant with tem- 
perature), frequency renormalization and lifetime effects, 
as is visualized in Fig. [5] LOPT predicts 



Po uo 



coth (OH) 
\2TJ 



suj ) 



s=± 



, (20) 



which breaks down completely at T ^ (uiq,u)q/T) even in 
the weak-coupling regime, leading to the complete sup- 
pression and even a sign change in G at high tempera- 
tures. 
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FIG. 2: (Color online) The electric conductance from 
bosonization (blue solid) and from LOPT (red dashed line) 
is shown in the left panel as a function of the temperature 
for weak electron- vibration coupling T — 0.2F2, W = Wioo =>■ 
uj±/ljo = ±0.89— iO. 02, signaling the limitations of the LOPT. 
The deviations grow with the coupling. The right panel shows 
the real and imaginary part of the Green's function together 
with the LOPT result at T = 0. 



Discussion. Experimentally, the renormalized vibra- 
tion frequency is determined from the current-voltage 
characteristics: the dl/dV exhibits a step down in the 
large contact broadening limit [17j at eV =Rew p ±, which 
sets the characteristic temperature range of the predicted 
conductance change. As a very rough estimate, the T = 
conductance is ~ pf (u> — * eV) , showing the multiphonon 
structures with the Franck Condon suppression [20j in the 
inset of Figs. [1] and [2 We mention, that low transmission 
junctions exhibit a step up in both the out-of-equilibrium 
dl/dV and the equilibrium spectral function at the exci- 
tation of the first vibrational mode. Temperature depen- 
dent transport is feasible on molecular devices [3l|, albeit 
large temperature variations are often accompanied by 
structural changes or the mechanical deformation of the 
junction. Therefore, the basic ingredients for the obser- 
vation of the predicted conductance step are the extreme 
mechanical stability of the device and a low enough fre- 
quency of the vibrational mode. Fullerenes like Cgo were 
found[6|] to possess a center-of-mass oscillation of 50 K, 
which can be lowered by considering heavier members of 
their family (e.g. Cuo). In this case, the oscillation fre- 
quency can accurately be estimated based on the interac- 
tions (electrostatic, van der Waals) between the molecule 
and the electrodes. In addition, the intercage vibrational 
modes of C140 start from 25 K0]. Another promising 
configuration, using a suspended quantum dot phonon 
cavity [32]| . possesses a vibrational frequency of 0.8 K. 

Although the electron- vibration coupling (g v ) is hardly 
controllable, the parameters A and po can be tuned by 
varying the contact DOS, which can drive the system 
towards stronger effective couplings. This decreases the 
temperature window for the conductance step. Given a 
low vibrational frequency, condition A 3> (T, o>q) is eas- 
ily met, therefore we expect that a dedicated setup al- 
lows the observation of the universal conductance steps, 
while low temperature dl/dV measurements can reveal 
the multiphonon structures in the local spectral func- 
tion on the molecule. Even for molecular devices, where 



the direct temperature dependence cannot be traced, our 
results demonstrate the role of multiphonon scattering 
processes in the room temperature conductance of the 
junction. 

In summary, we studied non-perturbatively the elec- 
tric and heat conductance through a single level coupled 
to a vibrational mode by mapping it to the Yu- Anderson 
model. With increasing temperature, the conductances 
drop to half of their T = value due to increasing inelas- 
tic scattering in the large contact broadening limit. We 
argue that its experimental observation is within reach 
in stable contacts with a low vibrational mode. 

This work was supported by the Hungarian Scien- 
tific Research Funds under grant numbers K72613 and 
K76010, and by the Bolyai program of the Hungarian 
Academy of Sciences. 
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